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A three-dimensional strong-topological-insulator or -semimetal hosts topological surface states 
which are often said to be gapless so long as time-reversal symmetry is preserved. This narrative can 
be mistaken when surface state degeneracies occur away from time-reversal-invariant momenta. The 
mirror-invariance of the system then becomes essential in protecting the existence of a surface Fermi 
surface. Here we show that such a case exists in the strong-topological-semimetal Bi 4 Se 3 . Angle- 
resolved photoemission spectroscopy and ab initio calculations reveal partial gapping of surface 
bands on the B^Ses-termination of Bi 4 Se 3 (lll), where an 85 meV gap along TK closes to zero 
toward the mirror-invariant FM azimuth. The gap opening is attributed to an interband spin-orbit 
interaction that mixes states of opposite spin-helicity. 
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Topological insulator materials (TIM) [1], which in¬ 
clude insulators [2-4] and semimetals [5-10], possess an 
inverted bulk band gap that hosts unusually robust, 
spin-helical topological surface states (TSS) at the ma¬ 
terial’s boundaries. The presence of TSS is guaran¬ 
teed by the topology of the bulk bands and the states 
can only be removed from the Fermi energy if protect¬ 
ing symmetries are broken. Strong topological insula¬ 
tor (STI) [11, 12] materials hold a special distinction, 
as they are often said to host gapless TSS protected by 
time-reversal symmetry (TRS) on every surface termina¬ 
tion of the crystal. We will show it is possible for the 
TSS of material in a STI phase to intrinsically acquire 
a finite gap at all momenta in the surface Brillouin zone 
(SBZ) which do not lie on a mirror-invariant azimuth. In 
that case, the mirror-symmetry (MS) of the crystal lat¬ 
tice must be intact to guarantee the existence of a sur¬ 
face Fermi surface. Earlier work by Teo, Fu, and Kane 
[13] had foreseen this possibility, however, no concrete 
examples have been obtained in experiment. Here we 
show that such a system exists in the strong-topological- 
semimetal BUSes. The finding of mirror-protected TSS 
in this system also provides direct confirmation that 
the strong-topological-insulator phase can coexist with 
the topological-crystalline-insulator phase of matter, as 
was suggested by Rauch et al [14] for the case of Bi- 
chalcogenides. 

BUSes is a rhombohedral superlattice material con¬ 
sisting of alternating Bi 2 and Bi 2 Se 3 layers stacked along 
the (111) direction. Previously, Dirac cone-like TSS were 
found in the gap between the first fully occupied bulk 
valence band (BVB) and the hole-like bulk conduction 
band (HBCB) on two different surface terminations of 
Bi 4 Se 3 (lll) [10]. It was determined that the TSS result 
from a parity inversion at the F-point of the bulk Bril¬ 
louin zone (BBZ), a characteristic shared with the Bi 2 Se 3 


parent compound [15]. Although effects of hybridization 
between surface and bulk electrons were mentioned, it 
was not made clear that the states of the upper Dirac 
cone appearing on the Bi 2 Se 3 -terminated surface, which 
have an electron-like dispersion, must have crossed the 
HBCB to reach the Fermi level. The possibility of new 
topological constraints on the surface electron-structure 
above the HBCB was also not explored. Earlier work [9] 
predicted that the band gap above the HBCB is char¬ 
acterized by a single parity inversion at the E-point of 
the BBZ, which demands that TSS within this gap must 
come in an odd number of pairs [13]. If that is the case, 
then an odd number of surface bands must appear within 
the gap to meet the state which crossed the HBCB. More¬ 
over, this parity inversion is away from the center of the 
BBZ and the crystal does not cleave at the center of in¬ 
version. Under these conditions, the TSS pairs are not 
constrained to have Dirac points at time-reversal invari¬ 
ant momenta (TRIM) [13]. 

Consistent with the prediction in ref. [9], it is in 
the momentum-space region outside the HBCB edge 
that evidence of partially gapped TSS is found through 
ab initio calculations and angle-resolved photoemission 
spectroscopy (ARPES) measurements on the Bi 2 Se 3 - 
terminated surface of Bi 4 Se 3 (lll). The TSS degenerate 
away from TRIM on the mirror-invariant FM azimuth 
in the SBZ. The bands become separated elsewhere, and 
an 85 meV gap between the surface state branches is 
measured along in ARPES. The origin of this gap is 
accounted for in a model for the spin-orbit interaction on 
a (111) crystal surface. These findings place 2D electron 
hybridization within the subject of pristine STI materials 
and provide direct evidence for MS protection of surface 
states in a Bi-chalcogenide. 

Single crystals of Bi 4 Se 3 were synthesized following a 
previously reported procedure [10]. ARPES was per- 
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formed using a Scienta SES -100 electron spectrometer 
at beamline 12.0.1 of the Advanced Light Source with 
a combined instrumental energy resolution of ^12 meV 
and an angular resolution better than ± O.OTA. The 
sample was cleaved under ultrahigh vacuum conditions 
(< 5.0 X 10“^ Pa) and kept at ^15 K. Temperature was 
measured using a silicon sensor mounted near the sample. 

Electron-structure calculations were performed in the 
framework of density functional theory (DET) using the 
WIEN2K code [16] with a full-potential linearized aug¬ 
mented plane-wave and local orbitals basis together with 
the Perdew-Burke-Ernzerhof [17] parametrization of the 
generalized gradient approximation, using a slab geom¬ 
etry. Experimentally determined lattice parameters and 
atom positions were used to construct the slabs. The 
plane-wave cutoff parameter RMxKniax was set to 7 and 
the Brillouin zone was sampled by 9 A; points. Spin- 
orbit coupling (SOC) was included. To study the Bi 2 Se 3 - 
terminated surface, a slab was constructed of 6 Bi 2 Se 3 
layers and 5 Bi 2 layers, with 10 A of vacuum between 
adjacent slabs. The contribution of the surface atoms to 
the overall surface electronic structure was determined by 
calculating the partial contribution of each atomic basis 
set to the wave functions at all k points. 

The calculated electron-structure for Bi 2 Se 3 - 
terminated slabs of BL 803 ( 111 ) is shown in Eigure 1 . A 
surface state crossing protected by MS only is indicated 
by the red circle and the blue circle indicates a crossing 
which is "dually protected" [14] by both MS and TRS. 
The TSS within the blue-circled region he between the 
BVB and HBCB. These TSS were the primary focus 
of previous investigations [9, 10]. The TSS that cross 
within the red-circled region, which are the focus of the 
present work, he above the HBCB. The parity invariants 
of the bulk band structure counted up to the HBCB 
were previously determined to be + 1 , + 1 , + 1 , and — 1 
at the T, Z, L, and E points of the BBZ, respectively 
[9]. The product of the parity invariants is —1, which 
characterizes the gap above the HBCB as a STI-type 
[ 11 , 12 ]. Applying the methods of ref. [13] to our case, an 
odd number of TSS pairs are expected to exist between 
the surface projections of E (M) and T (f). Indeed, 
we observe a single pair of TSS that cross each other 
along the TM azimuth and degenerate with different 
ends of the bulk gap at T and M. Group-theoretical 
considerations indicate why this crossing is allowed even 
while the surface states are seen to be gapped along the 
VK azimuth. 

At the ( 111 ) surface, the R3m symmetry of the crys¬ 
tal reduces to C 3 V. Eor wave-vectors lying between T 
and M, the point-group symmetry reduces to Cg, which 
contains two irreducible representations characterized by 
mirror eigenvalues of Through the definition of the 
mirror operation [13], it is easily shown that the two irre¬ 
ducible representations correspond to states of opposite 
spin-helicity, which cannot hybridize with each other on 
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Figure 1. Calculated band structure for_ the Bi 2 Sea- 
terminated slabs of Bi 4 Se 3 plotted along the MVK path in the 
surface Brillouin zone. The size of the circular plotting mark¬ 
ers indicate the contribution from the surface layer. Shaded 
regions indicate the projection of bulk electrons. The red cir¬ 
cle contains the mirror-symmetry-protected crossing of TSS 
and the red double arrow indicates the corresponding TSS 
avoided crossing gap A. The blue circle contains a TSS cross¬ 
ing protected by both TRS and MS. 


the mirror-invariant T M azimuth. This explains why the 
crossing circled in red is allowed and, indeed, protected 
by the crystal’s mirror symmetry. In contrast, the point 
group of the wave-vectors along VK is Ci. By symmetry, 
crossings between V and K are avoided, even for states 
of opposite spin-helicity. The same is true for all wave- 
vectors in the SBZ which do not he on a mirror-invariant 
azimuth. The double-arrow in panel (a) indicates what 
can therefore be understood as a hybridization gap A re¬ 
sulting from the avoided crossing of spin-helical surface 
states, as will be discussed later. Note that the combi¬ 
nation of C 3 and time-reversal symmetry imply that the 
surface Eermi surface will consist of six equivalent pock¬ 
ets that each enclose a surface state degeneracy point. 
This observation is consistent with the definition of a 
STI material as put forward by Teo, Eu, and Kane [13]. 

The "^-dependence of the gapped structure results 
from competing SOC interactions which couple to dif¬ 
ferent components of the spin degree of freedom. This 
can be captured in a model for the SOC Hamiltonian 
Hsoc ^ X using k • p theory. Polar coor¬ 

dinates are chosen with T as the origin and 0 denoting 
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the in-plane azimuthal angle from the /c^^-axis, aligned 
to the direction. If the energy-splitting of the TSS 
near k = 0 is taken to equal 2vhv^ then Hgoc to first- 
order in k couples the in-plane, tangential component of 
spin (at) to the out-of-plane electrostatic potential gra¬ 
dient with a strength v as Hi = v{k — K)(Jt. To 

third-order in A;, there appears a second term [18] that 
couples the out-of-plane component of spin to the in¬ 
plane crystalline potential gradient with a strength A as 
H 2 — ^k^cos{30)az. Together, 

Xk^cos{30) iv{k — 

—iv{k — —Xk^cos{30) 

and we find the spin-orbit contribution to the dispersion 

EsocCi^)± = 



where we refer to + and — as the upper Dirac branch 
(UDB) and lower Dirac Branch (LDB), respectively. At 
k = K, the magnitude of the gap between the DBs is de¬ 
termined solely by the second term under the square root. 
The gap carries the sign of an /-wave, which changes at 
the mirror-invariant azimuths, signifying a change in the 
^^-polarizations of the DBs (this could alternately be de¬ 
scribed as a crossing of bands with positive and negative 
z-polarization). The sign appears in the interband matrix 
element 


A = 


H2 (t) 


—2Xk^cos{30) 


for the spin-helical states typically invoked in the discus¬ 
sion of simple, gapless TSS. In this sense, the crystalline 
anisotropy should cause an interband SOC effect that 
gives rise to the gapped structure. The model only dif¬ 
fers from the TSS of the STI Bi 2 Te 3 [18] in the location 
of the Dirac point. No extra bands need to be inferred 
to achieve the complex, partially gapped TSS described 
below. 

Fig. 2 displays electron-structure calculated from the 
model Hamiltonian, for which we have chosen the val¬ 
ues V = 5.5 eVA, A = 55 eVA^, k. = 0.25 A'^ in 
rough approximation of the Bi 2 Se 3 -termination electron- 
structure. The reader should note the many omissions 
from this model such as spin-orbital entanglement (which 
will disallow a pure spin-eigenstate character for the TSS 
and limit the magnitude of spin-polarization) [19], and 
higher-order interaction terms [20]. Fig. 2(a) shows the 
spin-resolved band structure along the 0 = 0 azimuth, 
with the bands corresponding to the case A = 0 plotted 
in black. The minimum energy gap 6 is located inside of 
k = K, and we find it is the case that dit ks < 1 ^ for all 6>, 
as shown in panel (c). It is telling to inspect the in-plane 
helical spin-polarization, indicated by color scale, rela¬ 
tive to the magnitude of the out-of-plane polarization, 
indicated by marker size. For A; < 0.1 A'^, the spin is he¬ 
lical, but as the contribution of H 2 grows, the deviation 




Figure 2. Electron-structure given by the model Hamiltonian 
for parameters stated in the text: Band structure near the 
anti-crossing point along FK {ky = 0) direction (a). The color 
scale (a, inset) indicates the spin-polarization in y direction, 
while marker size indicates the degree of z-spin-polarization. 
Constant energy contours (CECs) of the upper Dirac branch 
(b). The color scale (b, inset) indicates the energy E of 
each contour measured with respect to the Dirac point. The 
momentum-space contour of the gap minimum overlayed onto 
the CEC at the Lifshitz transition energy is shown in (c). 


from linear dispersion is accompanied by increasing 2 ;- 
polarization. At ks the contribution of H 2 overtakes that 
oi Hi. AF k = the 2 ;-polarization is 100% and then at¬ 
tenuates to ^90% for A: > /^, where the spin-helicity has 
reversed. The scenario of competing SOC interactions 
takes place throughout the SBZ, resulting in unusual 
constant energy contour (CEC) shapes and topologies, 
displayed in panel (b). If the electron-dispersion were 
determined solely by the SOC in this model, there would 
exist two Lifshitz points [21] in the chemical potential lo¬ 
cated at /i = ±E{ks,0 = 0). At these points, the Fermi 
surface topology changes from six pockets enclosing the 
TSS degeneracies to one electron-pocket and one hole- 
pocket enclosing F. Near a Lifshitz point, straight edges 
in the CECs are centered on the FM direction, opposite 
of what would be expected for simple, gapless TSS on a 
(111) surface [18]. This same pattern appears in CECs 
probed by ARPES, described in Fig.3(c) below. 

Fig.3 shows the ARPES spectra of the terraced sur¬ 
face of Bi 4 Se 3 (lll) collected using 70 eV photons. Pre¬ 
vious photoemission electron microscopy studies revealed 
that surface is terminated by Bi 2 and Bi 2 Se 3 layers which 
are present in approximately equal proportion [10]. The 
Fermi surface in panel (a) has three distinct branches 
enclosing F. It was previously determined that the cir¬ 
cular, innermost branch is derived from bulk conduction 
electrons and Bi 2 -surface electrons, while the hexagonal 
branch outside of that is derived from Bi 2 -surface elec- 
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Figure 3. ARPES electron-structure near the center of the surface Brillouin zone: The Fermi surface (a,b), and constant 
energy contour at -185 meV (c). Band structure parallel to FM direction for = 0 (d) and kx — 0.21 A“^(i). Band structure 
parallel to f K direction at ky = —0.255 A“^ (e), ky = 0.255 A“^ (f), and ky = 0 (h). Energy distribution curves at the degeneracy 
point {kx, ky) = (0, —0.255) A“^ (g) and the saddle point of the lower Dirac branch {kx, ky) = (0.21, 0) A“^ (j). The momentum- 
space locations for panels (d-e) and (h-i) are indicated by lines overlayed onto the Fermi surface in (a). The gray-scale varies 
from white-to-black following the minimum-to-maximum photoemission intensity within each image individually. 


trons [10]. The outermost branch is derived from sur¬ 
face electrons of the Bi 2 Se 3 -termination, which display 
a 3-fold enhancement of intensity due to ARPES matrix 
element effects. The degeneration of TSS on the mirror- 
invariant FM {kx = 0) line in the vicinity of ky = ±0.255 
A'^ at Eb = —0.05 eV is observed in the band structure 
shown in panels (d-f). Crucially, the band structure at 
ky = ±0.255 A'^ (e) and ky = —0.255 A'^(f) appears 
to be identical, confirming that these bands consist of 
spin-polarized surface states, which have a 6-fold sym¬ 
metry on the underlying 3-fold-symmetric lattice as re¬ 
quired by TRS. Comparing the Fermi surface in panel 
(b) to that of the -185 meV CEC (Fig.3(c)), we observe 
the formation of a "teardrop-shaped" CEC that possess 
a straight edge centered on FM, similar to what is pre¬ 
dicted by the model Hamiltonian. Panels (h-j) reveal a 
saddle-point in the LDB near {kx,ky) = (0.21, 0) A'^ (in¬ 
dicated by red dashed-lines), where the minimum energy 
separation 6 between the DBs reaches a maximum value 
(with respect to the in-plane azimuthal angle) of 85 meV. 
Comparing (h-i), it is clear that the LDB reaches a min¬ 
imum with respect to kx and a maximum with respect to 
ky at this point. This would yield a Van Hove singular¬ 
ity [22] in the density of states. Interestingly, the UDB 
is at a local minimum with respect to both variables at 
the same point in momentum-space. The model Hamil¬ 
tonian, which retains particle-hole symmetry, predicts a 
saddle point in both DBs. 

The observation of gapped surface states on a (111) 
surface with large SOC is not without precedent, having 
been previously identified in heterostructures with BiAg 2 
surface alloys [23], however, the mechanism for the "in¬ 
terband SOC" between the antiparallel, spin-helical sur¬ 
face states in that case was left unspecified. The two- 
term model Hamiltonian approach shown in this work 


could be extended to describe not only BiAg 2 surface 
states, which are known to have a sizable coupling to the 
in-plane crystal potential gradient [24], but also many 
other systems, whether topologically trivial or not, in 
which spin-helical states intersect away from Kramer’s 
momenta, such as Pb quantum wells [25]. TIM, rather 
than topologically trivial materials, may offer more ro¬ 
bust platforms for studying this type of spin-gap physics 
in 2D electron systems. What is lacking at this time is a 
straightforward and reliable way of predicting if a given 
STI surface will possess gapped TSS. The present results 
indicate a need to merge different conceptualizations of 
the topological insulator (as defined in the abstract sense 
band topology [26]) in order to determine the conditions 
necessary for this phenomenon. 

Some have pointed out [13, 14, 27] that several well- 
known TIM possess bulk electron structure that can 
be characterized as topologically non-trivial using sep¬ 
arate methods. The parity invariant method is used to 
characterize Z 2 topological insulators [11, 12], which in¬ 
clude STIs, while mirror topological crystalline insulators 
(TCIs) are characterized by a non-zero difference in the 
number of counterpropogating "edge states" correspond¬ 
ing to a crystallographic mirror plane [13, 28, 29]. Teo, 
Fu, and Kane [13] had considered that the crossing of TSS 
at non-TRIM was a possibility for a strong Z 2 topological 
insulator surface, which motivated them to develop the 
foundational theory for TCIs. This Letter has presented 
an experimentally realized case in which the concepts 
of Z 2 topological insulators and TCIs have become en¬ 
tangled beyond precedent; completely breaking the MS 
would allow the TSS to become fully gapped, even while 
TRS remains unbroken. Surely, the significance of the Z 2 
topology in determining the electronic physics at all of 
the possible surfaces of a STI should be revisited. 
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